1/1 


AD-A138  752 
UNCLASSIFIED 


USAF  -  SCEEE  (SOUTHEASTERN  CENTER  FOR  ELECTRICAL 
ENGINEERING  EDUCATION)  S..(U)  SOUTHEASTERN  CENTER  FOR 
ELECTRICAL  ENGINEERING  EDUCATION  INC  S. .  AG  TURNER 
15  AUG  80  FUSRL-JR-80-0004  F49620-79-C-0038  F/G  19/1  NL 


microcopy  resolution  test  chart 

NATIONAL  BUREAU  Of  STANDARDS- 1963  4 


ADA1  38752 


1980  USAF  -  SCEEE  SUMMER  FACULTY  RESEARCH  PROGRAM 
Sponsored  by  the 

AIR  FORCE  OFFICE  OF  SCIENTIFIC  RESEARCH 
Conducted  by  the 

SOUTHEASTERN  CENTER  FOR  ELECTRICAL  ENGINEERING  EDUCATION 

FINAL  REPORT 

A  MODEL  FOR  THE  THERMAL  DECOMPOSITION  OF  TNT; 
THEORETICAL  REACTION  PROFILES 


Prepared  by: 
Academic  Rank: 


Professor  Almon  G.  Turner 
Professor 


fj  5/f  L  '  -WW 


Department  and  Department  of  Chemistry 

University:  University  of  Detroit 

Research  Location:  Frank  J.  Seiler  Research  Laboratory, 

USAF  Academy,  Colorado 


USAF  Research 
Colleague : 

Date: 

Contract  No: 


Dr.  Larry  P.  Davis 
August  15,  1980 
F49620-79-C-0038 


DT:r 

;  ... 

MAR  7  1984 

o 

LU 

— 1 

u_ 

This  document  »«..  L-en  appiovsd 

g 

j  1.1  public  reltasc  and  sale;  its 
!  oi  :.-jbution  unlin'/cd.  _ 

8  4  o  Z  \  o 

Acknowledgements 


The  author  would  like  to  thank  the  Air  Force  Systems  Command,  the 
Air  Force  Office  of  Scientific  Research  and  the  Southeastern  Center  for 
Electrical  Engineering  Education  for  financial  support  via  a  Summer 
Faculty  Research  Program  Fellowship.  He  should  like  to  acknowledge  the 
hospitality  extended  to  him  by  the  members  of  the  Frank  J.  Seiler  Research 
Laboratory,  USAF  Academy. 

He  should  particularly  like  to  acknowledge  the  cooperation  of  Capt 
Larry  P.  Davis  for  introducing  him  to  the  area  of  the  chemical  decomposi¬ 
tion  of  explosives  and  for  the  many  helpful  discussions  held  in  this  area. 


A  MODEL  FOR  THE  THERMAL 


DECOMPOSITION  OF  TNT;  THEORETICAL 
REACTION  PROFILES 
by 

Almon  G.  Turner 

/ 

ABSTRACT 


^The  molecule  1-nitropropylene  has  been  investigated  as  a  model  system 
for  the  simulation  of  the  thermal  decomposition  of  TNT.  Two  distinct  types 
of  reaction  mechanisms  were  considered:  Intramolecular  Mechanisms  and 
Bimolecular  Mechanisms.  Intramolecular  mechanisms  investigated  include 
an  oxidative  hydrogen  atom  transfer  to  form  the  diradical  aci  nitropropylene 
and  an  oxidative  oxygen  atom  insertion  reaction  to  form  1-nitro- 3- hydroxy 
propylene.  Semi  empirical  molecular  orbital  calculations  (MNDO  and  MINDO/3) 
were  carried  out  to  obtain  a  reaction  profile  for  these  mechanisms,  and 
indicated  that  the  oxygen  atom  insertion  reaction  should  lead  to  the 
reaction  products,  methyl  nitrite  and  acetylene.  This  is  not  in  accord 
with  experiment.  The  profile  for  the  oxidative  hydrogen  atom  transfer  was 
found  to  reproduce  many  of  the  features  known  for  the  thermal  decomposition 
of  TNT. 

Bimolecular  mechanisms  considered  included  an  oxidative  hydrogen  atom 
transfer  from  one  nitropropylene  molecule  to  another  to  form  the  1-nitro¬ 
propylene  radical  and  the  aci  form  of  nitropropylene  and  an  intermolecular 
oxidative  insertion  reaction  to  form  l-nitro-3-hydroxy  propylene  and  1- 
nitroso  propylene.  Preliminary  reports  are  given  for  these  bimolecular 
mechanisms. 


I.  INTRODUCTION: 

A  large  number  of  experimental  studies  have  been  carried  out  in  an 
effort  to  obtain  an  understanding  of  the  thermal  decomposition  of  trinitro¬ 
toluene  (TNT) These  studies  have  resulted  in  the  proposal  of  four  dif- 

(2) 

ferent  possible  reaction  mechanisms  for  the  process  .  The  mechanisms 
proposed  are: 

a.  Intramolecular  Oxidative  Hydrogen  Atom  Transfer. 


X  X 


(I)  (T I) 

b.  Intramolecular  Oxidative  Oxygen  Atom  Transfer. 


(Ill) 


V 


Recent  advances  in  semiempirical  molecular  orbital  methods  (MNDO,  MINDO/3) 

O 

have  demonstrated  that  the  standard  enthalpies  of  formation,  AH  ,  can  be 

(3)  f 

calculated  with  some  confidence  for  small  molecules  ,  These  same 

methods  have  been  utilized  to  study  reaction  mechanisms  for  elementary 

(4) 

processes  of  the  type  which  occur  in  the  above  mechanisms  .  The  shear 
size  of  the  TNT  molecules  prohibits  an  exhaustive  investigation  of  the 
proposed  mechanisms  by  the  MNDO  method.  Accordingly  it  would  be  advanta¬ 
geous  to  find  a  smaller  molecular  system  which  might  possess  the  salient 
features  of  TNT,  and  could  thus  be  used  to  model  the  mechanisms  proposed 
above.  The  author  has  suggested  that  the  1-nitropropylene  molecule  might 
be  suitable  for  these  purposes.  The  research  reported  below  is  concerned 
with  the  use  of  1-nitropropylene  to  model  the  thermal  decomposition  of  TNT. 

II.  OBJECTIVES 

The  main  purpose  of  this  research  was  to  ascertain  the  extent  to 
which  mechanisms  analogous  to  those  listed  above  but  studied  in  the 
system  1-nitropropylene  can  be  used  to  model  the  decomposition  of  TNT. 
Specific  objectives  were: 

■i  (1)  To  study  the  thermodynamics  of  the  proposed  reactions  and  compare 

the  results  to  the  corresponding  values  for  the  TNT  system. 

(2)  To  obtain  reaction  profiles  for  each  of  the  proposed  mechanisms  and 

to  identify  the  transition  state  for  each  of  the  mechanisms.  Success¬ 
ful  identification  of  a  transition  state  enables  one  to  calculate 

f 

1 


the  energy  of  activation  for  the  process,  which  in  turn  can  be  compared 
to  results  obtained  from  the  experimental  study  of  the  decomposition. 


III.  THEORETICAL  METHODS 

The  principal  method  employed  to  study  the  decomposition  of  1-nitropro- 

pylene  was  the  MNDO  (Modified  Neglect  of  Differential  Overlap)  molecular 

(34) 

orbital  method  of  Dewar  et.  al.  ’  .  The  method  was  applied  in  its  "closed 
shell"  form  to  those  mechanisms  which  did  not  involve  radical  species,  while 
it  was  combined  with  the  method  of  configuration  interaction  to  treat  radi¬ 
cal  species^.  3x3  configuration  interaction  was  carried  out  using  the 
ground  state  configuration  and  the  two  configurations  obtainable  from  it 
using  the  two  lowest  unoccupied  molecular  orbitals  obtained  from  the  ground 
state  calculation. 

IV.  THERMODYNAMIC  STUDIES 

The  specified  chemical  reaction  (mechanisms)  studied  were: 

i.  Intramolecular  Hydrogen  Atom  Transfer 

ch3chchno2  - »  -ch2chchno2h 

1-nitropropylene  (I)  ACI-nitropropylene  (II) 

ii.  Intramolecular  Oxygen  Atom  Transfer 

ch3chchno2  - ►  hoch2chchno 

l-nitroso-3-hydroxy  propylene  (III) 

ch3chchno2h 

(V) 

1-nitrosopropylene  (VI) 

+  hoch2chchno2 

l-nitro-3-hydroxy  propylene  (VII) 

The  ground  state  energies  and  heats  of  formation  were  calculated  for 
species  I  -  VII  and  the  corresponding  heats  of  reaction  calculated.  The 
1-nitropropylene  molecule  was  calculated  in  a  geometry  where  the  nitro 


iii.  Intermolecular  Hydrogen  Atom  Transfer 

2  ch3chchno2  - -ch2chchno2  + 

(IV) 

iv.  Intermolecular  Oxygen  Atom  Transfer 

2  ch3chchno2  - ►  ch3chchno 


group  was  rotated  about  the  carbon  atom  so  as  to  have  it  in  the  same  con¬ 
formation  that  it  would  have  were  it  a  part  of  a  TNT  molecule.  The  energy 
of  the  true  ground  state  form  (unrotated)  was  found  to  be  18.10  Kcal/mole. 
The  energy  of  the  rotated  "TNT  like"  form  was  calculated  to  be  18.99  Real/ 
mole.  The  results  of  these  calculations  are  shown  below. 

Calculated  Standard  Heats  of  Formation 


Mechanism 

Reactant 

Product(s) 

AHf 

TNT  Case 

i 

18.99 

33.50 

+14.51 

16.6 

ii 

18.99 

-30.32 

-49.31 

-49.9 

ill 

37.98 

IV  48.63;  V  8.88 

19.53 

16.8 

iv 

37.98 

VI  12.45;  VII  022.98 

-48.51 

-52.4 

The  experimental  activation  energy  for  the  initiation  step  in  the 

O 

decomposition  of  TNT  is  known  to  be  40-46  Kcal/mole.  Since  the  AH^  for  the 
reaction  is  a  lower  bound  for  the  corresponding  activation  energy  all  the 
above  mechanisms  are  possible.  The  corresponding  enthalpy  changes  have  been 
calculated  for  three  of  the  four  mechanisms  in  the  case  of  TNT.  The  values 
are  given  in  the  last  column  of  the  table.  Note  the  relative  order  of 
agreement  is  good.  On  the  basis  of  thermodynamics  the  1-nitropropylene 
system  seems  to  reflect  the  same  energetic  behavior  as  the  TNT  system. 

V.  INTRAMOLECULAR  HYDROGEN  ATOM  TRANSFER  MECHANISM 

Optimization  of  the  ground  state  energy  of  aci  nitropropylene  yielded 
a  value  of  33.50  Kcal/mole  for  the  heat  of  formation  of  the  product  molecule. 
The  hydrogen  atom  transfer  was  initially  viewed  by  using  the  hydrogen  atom 
oxygen  atom  distance  as  the  reaction  coordinate.  The  distance  varies  from 

O  O 

3.73  A  in  1-nitropropylene  to  0.95  A  in  product  aci  form.  A  reaction  pro¬ 
file  was  calculated  along  this  path  at  values  for  the  reaction  coordinate 
of  3.73,  3.27,  2.77,  2.27,  1.77,  1.50,  1.40,  1.35,  1.30,  1.25,  1.20,  1.175, 

O 

1.15,  1.125,  1.10,  1.00,  0.95  A.  A  maximum  was  obtained  in  the  vicinity 

O 

dQ_H  =  1.20  where  AH^  =  81.4  Kcal/mole,  which  corresponds  to  an  activation 
energy  of  62.4  Kcal/mole.  The  product  state  was  identified  at  a  heat  of 
formation  of  about  27.4  Kcal/mole  and  did  not  correspond  geometrically 


to  the  aci  form  of  1-nitropropylene.  Attempts  were  made  to  start  from  this 
new  product  state  and  back  calculate  to  the  initial  state.  These  attempts 
were  unsuccessful. 


The  next  series  of  calculations  performed  consisted  of  starting  the 
reaction  profile  in  the  product  state  and  attempting  to  find  a  path  which 
would  lead  to  the  starting  material.  This  was  done  by  defining  the  reac¬ 
tion  coordinate  to  be  the  hydrogen  atom  -  methyl  carbon  atom  distance. 

O  O 

This  distance  varies  from  3.234  A, to  1.497  A  in  the  reactant.  Points 
were  calculated  at  values  of  3.20,  3.00,  2.50,  2.00,  1.75,  1.60,  1.50,  1.40, 

O 

1.35,  and  1.30  A.  A  profile  was  obtained  with  a  maximum  energy  of  86.3 
Kcal/mole  corresponding  to  an  activation  energy  of  52.8  Kcal/mole  occurring 

o 

at  d  =  1.37  A.  An  attempt  was  made  to  reverse  the  calculation  and  it 
was  found  that  the  path  was  not  reversible,  the  reverse  path  leading  to  a 
new  maxima,  in  energy  in  excess  of  76  Kcal/mole.  The  two  profiles  inter- 

O 

sected  at  about  1.59  A  with  an  energy  of  about  69  Kcal/mole.  Application 
of  the  MINDO/3  method  led  to  the  same  basic  results.  From  these  calcu¬ 
lations  it  was  concluded  that  neither  the  hydrogen  -  oxygen  atom  distance 
nor  the  hydrogen  atom  -  methyl  carbon  atom  distance  was  a  suitable  reaction 
coordinate.  It  is  probably  necessary  to  use  both.  In  addition  as  the 
molecule  proceeds  from  reactant  to  product  a  considerable  readjustment 
occurs  in  the  H„C. -C  =C  -N .  framework.  The  C,-C0  and  C  -N.  bonds  are 
shortened  considerably,  while  a  lengthening  occurs  in  the  C2~C3  b°nd- 

Accordingly,  a  series  of  calculations  were  carried  out  for  d  distances 

U-H 

of  0.96,  1.00,  1.10,  1.20,  1.70,  1.40,  and  1.50  where  the  C^-C2=C^-N 
framework  was  fixed  at  values  which  corresponded  to  one  third,  one  half, 
and  two  thirds  of  the  way  along  the  progression  from  starting  material 
to  product  material.  The  results  of  these  calculations  are  in  Table  1. 

O 

Listed  are  the  AH^'s  and  the  algebraic  sign  of  dE/^ 

t 

i 

i 

f 

I 


Table  l  GRID  MOVING  FROM  REACTANT  TO  PRODUCT 


0-H 


!  0.96 

1.00 

1.10 

1.20 

1.30 

1.40 

1.50 

Starting 

Material 

r* 

62.6- 

63.8+  V 

\91.6- 

81.8- 

71.5- 

61.8- 

53.0- 

1/3 

42.9- 

43.7+ 

52?0>> 

83.1- 

73.4- 

64.0- 

55.4- 

1/2 

76.6- 

37.4+ 

45.7+ 

85.8- 

76.8- 

67.8- 

59.3- 

2/3 

32.6- 

33.7+ 

45.6+ 

89.3- 

81.2- 

72.4- 

64.0- 

Product 

Material 

27.6- 

28.4+ 

36.9+ 

„  96.8- 

88.9- 

81.1- 

The  line  delineates  the  region  of  the  hypersurface  where  maxima  are 
located. 


The  largest  maxima  being  in  the  region  of  120(d  Xl.30  with  the  geometry 

being  virtually  TNT  or  product  like.  Calculations  were  then  performed  to 

test  the  reversibility  of  the  reaction  coordinate  d  .  Starting  with  the 

TNT  like  geometry  (two  thirds)  at  d.  =  1.10,  d„  ,,  was  allowed  to  increase 

U— rt  O-H 

to  1.20,  1.30  and  1.40.  At  each  point  a  lower  energy  form  was  found.  This 

indicated  that  the  =  framework  must  also  be  relaxed. 

A  grid  search  was  conducted  by  performing  calculations  which  start 

with  a  TNT  like  geometry  and  allowing  d  to  vary  from  1.00  to  1.60  and 

0— H 

at  the  same  time  allowing  d„  u  to  vary  from  1.20-1.80.  The  results  are 

L— H 

shown  in  Table  2. 


Calculations  wore  next  carried  out  by  starting  with  tbe  product  like  geome¬ 
tries  for  d  equal  to  1.00  and  attempting  to  proceed  to  the  TNT  like 
geometries  for  each  C-H  distance.  It  appears  that  some  points  of  lower 
energy  can  be  located  in  this  manner. 

In  summary,  the  location  of  transition  state  for  this  mechanism  is 
particularly  ditficult  and  involves  a  great  deal  of  trial  and  error  calcu¬ 
lation.  It  is  felt  that  a  transition  state  is  being  approached  asymptoti¬ 
cally  but  the  calculation  points  to  the  severe  need  for  an  analytical 
approach  to  the  problem. 

VI .  INTRAMOLECULAR  OXYGEN  ATOM  TRANSFER  MECHANISM 

A  calculation  of  the  optimized  ground  state  energy  of  the  product 
molecule,  l-nitroso-3-hydroxy  propylene  yielded  a  standard  enthalpy  of 
formation  of  -30.32  Kcal/mole.  The  methyl  carbon-oxygen  atom  distance 

O 

was  employed  as  a  reaction  coordinate  and  it  was  allowed  to  vary  from  3.73  A 

G 

(starting  material)  to  1.40  A  (product  material).  A  reaction  profile  was 
obtained  which  showed  an  activation  energy  of  96  Kcal/mole  and  a  maximum 

O  O 

present  at  d^  ^  =  1.75  A.  An  examination  of  the  product  state  (AHj.  = 

21.71  Kcal/mole)  showed  it  to  consit  of  molecular  acetylene,  C^H^;  and 
methyl  nitrite,  CH^ONO.  A  rerun  of  the  calculation  using  configuration 
interaction  led  to  new  transition  state  of  a  higher  heat  of  formation  than 
that  calculated  previously.  As  such,  it  was  not  pui  jued  further.  An  attempt 
was  made  to  establish  a  reaction  profile  in  the  direction  product  to 
reactant.  It  also  led  to  a  higher  activation  energy. 

VII.  INTERMOLECULAR  MECHANISMS 

Work  was  just  being  initiated  on  these  mechanisms  as  the  period  of  my 
fellowship  came  to  end.  Planary  calculations  were  begun  on  a  mechanism  in 
which  the  1-nitropropylene  molecules  were  stacked  end  to  end  with  the  hydro¬ 
gen  atom  to  be  transferred  of  one  molecule  located  symmetrically  between 
the  two  oxygen  atoms  of  the  nitro  group  of  tb»  second  molecule.  The  distance 
from  the  oxygen  atom  of  the  second  molecule  to  the  methyl  carbon  atom  of  the 
first  molecule  was  used  as  a  reaction  coordinate.  These  efforts  are  continu¬ 
ing  at  the  present  time. 


VIII.  AN  ANALYTICAL  APPROACH 


The  principal  difficulties  encountered  in  these  studies  seems  to  be 
the  recognition  of  what  to  use  for  a  reaction  coordinate  in  a  given  mechan¬ 
ism.  The  reaction  coordinate  is  generally  not  a  simple  atom-atom  distance, 
bond  angle  or  dihedral  angle  but  rather  it  is  expressible  as  a  linear 
combination  of  those  atom-atom  distances,  bond  angles  and  dihedral  angles 
which  change  markedly  as  one  proceeds  from  reactant  to  product.  As  such, 
it  can  be  uniquely  identified  once  it  is  recognized  what  internal  coordinates 
(bond  distance  -  bond  angle  basis)  enter  into  it.  One  way  to  proceed  would  be 
to  do  a  few  planary  calculations  to  identify  which  internal  coordinates, 

| q |  enter  the  reaction  coordinate,  $  and  then  solve  the  required  set  of 
linear  equations  to  identify  the  reaction  coordinate.  Explicitly,  let 

$  =  E  qt  i  =  1,  n  (1) 


Consider  the  energy  as  a  function  of  <J> , 

E  =  E  (<f)  (2) 

Since  we  are  always  seeking  a  transition  state,  E  is  a  maximum  for  such 
a  state, we  can  write 

E  =  X<p2  +  Y$  +  1  (3) 

and  we  know  that  X  <0  and  Y  >  0.  The  reaction  coordinate  must  belong 
to  one  of  the  irreducible  representations  of  the  point  group  which  is 
appropriate  to  the  transition  state  geometry. 

2  i  2 

From  a  set  of  n  values  for  the  jqj  one  can  calculate  n  values  of  the 
energy  e.  Since  the  energy  is  a  quadratic  function  of  the  |q|  we  can  write 

‘  '  blnqin  +  b10qil  +  b20ql2  +  bnoqin  +  bl 


£1  "  bllqil2  +  b12ql22  + 


(4) 


'n2 


=  bllqnl 


2  2 

b-„q  ,  +  b,  q  +  b, _q  ,  +  b  q 

12nn2  In  n,n  10  n,l  no  n,n 


+  b. 


The  linear  equations  (4)  can  be  solved  to  obtain  the  coefficients  |b| . 
These  coefficients  in  turn  are  related  to  coefficients  | a  f of  equation  (1) . 
Namely 

for  quadratic  terms 


for  linear  terms 


b 


io 


Y 


a  i 


(6) 


and 

bj  =  Z  (7) 

From  eqn.  (6)  we  in  effect  know  the  coefficients  jaj  to  within  the 
normalization  constant  Y.  If  we  redefine  the  norm  of  the  space  spanned 
by  the  reaction  coordinate,  we  can  calculate  the  jaj  and  then  obtain  our 
reaction  coordinate. 


IX.  RECOMMENDATIONS 

On  the  basis  of  the  above  research  the  following  recommendations  can  be 

made: 

1.  That  the  grid  search  method  employed  in  part  V  be  continued  until  a 
transition  state  is  obtained  which  will  reversibly  interconnect  reaction 
and  product  for  the  intramolecular  hydrogen  atom  transfer  mechanism. 

2.  That  the  intramolecular  oxygen  atom  transfer  be  rejected  as  a  possible 
reaction  path  in  the  initiation  of  the  thermal  decomposition  of  TNT. 

3.  That  the  analytical  approach  (part  VIII)  to  finding  a  reaction 
coordinate  be  tested  in  the  context  of  a  known  reaction  mechanism  and  if 
successful  then  used  to  pursue  the  study  of  the  bimolecular  mechanisms  of 
part  VII. 

4.  That  the  questions  concerned  with  the  effects  of  overcompleteness  and 
undercorapleteness  in  the  choice  of  the  set  of  variables  |q[  be  investi¬ 
gated  on  the  analytical  approach  (part  VIII). 
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